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With the experimental data from STAR and PHENIX on the centrahty dependence of the pr 
spectra of protons and anti-protons produced at mid-rapidity in Au-|-Au coUisions at 200 GeV, 
we show that for protons and anti-protons there exists a scahng distribution independent of the 
colhding centrahty. The scahng functions can also describe data from BRAHMS for both proton 
and anti-proton spectra at j/ = 2.2 and 3.2. The scaling behaviors are shown to be incompatible 
with the usual string fragmentation scenario for particle production. 
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I. INTRODUCTION 



One of the most important quantities in investigating 
properties of the medium produced in high energy col- 
lisions is the particle distribution for different species of 
final state particles. RHIC experiments have found a lot 
of novel phenomena from the particle spectra, such as the 
unexpectedly large p/n ratio at pt ~ 3 GeV/c the 
constituent quark number scaling of the elliptic flows , 
and strong nuclear suppression of the pion spectrum in 
central Au-f Au collisions etc. From the spectrum one 
can learn a lot on the dynamics for particle production. 

In many studies, searching for a scaling behavior of 
some quantities vs suitable variables is useful for unveil- 
ing potential universal dynamics. A typical example is 
the proposal of the parton model from the x-scaling of 
the structure functions in deep-inelastic scatterings 
Quite recently, a scaling behavior [S] of the pion spectrum 
at mid-rapidity in Au-|-Au collisions at RHIC was found, 
which related spectra with different collision centralities. 
In Q the scaling behavior was extended to non-central 
region, up to 77 = 3.2 for both Au-|-Au and d-|-Au colli- 
sions. The same scaling function can be used to describe 
pion spectra for pt up to a few GeV/c from different col- 
liding systems at different rapidities and centralities. The 
shape of pion spectrum in those collisions is determined 
by only one parameter (pt), the mean transverse momen- 
tum of the particle. It is very interesting to ask whether 
similar scaling behaviors can be found for spectra of other 
particles produced in Au+Au collisions at RHIC. In this 
paper, the scaling property of the spectra for protons and 
anti-protons is investigated and compared with that for 
pions. 

The organization of this paper is as follows. In Sec. II 
we will address the procedures for searching the scaling 
behaviors. Then in Sec. Ill the scaling properties of the 
spectra for protons and anti-protons produced in Au-f Au 
collisions at RHIC at ^ysNN — 200 GeV will be studied. 
We discuss mainly the centrahty scaling of the spectra at 
mid-rapidity and extend the discussion to very forward 
region with rapidity y — 2.2 and 3.2 briefly. Sec. IV is for 
discussions on the relation between the scaling behaviors 
and the string fragmentation scenario. 



II. METHOD FOR SEARCHING THE SCALING 
BEHAVIOR OF THE SPECTRUM 

As done in , the scaling behavior of a set of spectra 
at different centralities can be searched in a few steps. 
First, we define a scaled variable 



z = pt/K , 



and the scaled spectrum 
$(z) = A 



2TTpTdpTdy 



(1) 



(2) 



Pt—Kz 



with K and A free parameters. As a convention, we 
choose K — A — 1 for the most central collisions. With 
this choice $(z) is nothing but the pt distribution for 
the most central collisions. For the spectra with other 
centralities, we try to coalesce all data points to one curve 
by choosing proper parameters A and K. If this can 
be achieved, a scaling behavior is found. The detailed 
expression of the scaling function depends, of course, on 
the choice of A and K for the most central collisions. 
This arbitrary can be overcome by introducing another 
scaling variable 



u = z/ (z) = pt/ {pt) , 
and the normalized scaling function 

/>oo 

^(u) = / <S>{z)zdz 



Here (z) is defined as 

(z) = / z^{z)zdz/ I ^{z)zdz 



(3) 
(4) 

(5) 

By definition, 'i'{u)udu =^ u'i'{u)udu = 1. This 
scaled transverse momentum distribution is in essence 
similar to the KNO-scaling 7] on multiplicity distribu- 
tion. 



III. SCALING BEHAVIORS OF PROTON AND 
ANTI-PROTON DISTRIBUTIONS 

Now we focus on the spectra of protons and anti- 
protons produced at mid-rapidity in Au-|-Au collisions at 
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^/sWn = 200 GeV. STAR and PHENIX Collaborations 
at RHIC published spectra for protons and anti-protons 
at mid-rapidity for a set of colliding centralities [1, 0. 
STAR data have a.pT coverage larger than PHENIX ones. 
As shown in Fig. [U all data points for proton spectra at 
different centralities can be put to the same curve with 
suitably chosen A and K, by the procedure explained 
in last section. The parameters are shown in Table I. 
Except a few points for very peripheral collisions (cen- 
tralities 60-92% for PHENIX data and 60-80% for STAR 
data), all points agree well with the curve in about six 
orders of magnitude. The larger deviation of data at 
centralities 60-92% for PHENIX and 60-80% for STAR 
from the scaling curve may be due to the larger central- 
ity coverage, because the size of colliding system changes 
dramatically in those centrality bins. For simplicity we 
define v — ln(l -I- z), and the curve can be parameterized 
as 

$p(z) = 0.052 exp(14.9w - 16. 21;^ + 3.3i;^) . (6) 
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FIG. 2: Scaling behavior of the spectrum for anti-protons 
produced at mid-rapidity in Au-|-Au collisions at RHIC. The 
data are taken from 8', 9] . Feed-down effects are not corrected 
in the STAR data for p. The solid curve is from Eq. (7). 
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FIG. 1: Scaling behavior of the spectrum for protons pro- 
duced at mid-rapidity in Au-f-Au collisions at RHIC. The data 
are taken from [3, [y]. Feed-down corrections are considered 
in the data. The solid curve is from Eq. (6). 

Similarly, one can put all data points for anti-proton 
spectra at different centralities to a curve with other sets 
of parameters A and K which are given also in TABLE I. 
The agreement is good, as can be seen from Fig. [21 with 
only a few points in small px region for peripheral colli- 
sions departing a little from the curve. For anti-proton 
the scaling function is 

$p(z) = 0.16 exp(13w - + 2.9v^) , (7) 

with V defined above. 

To see how good is the agreement between the fitted 
curves in Figs. 1 and 2 and the experimental data, one 
can calculate a ratio 

B = experimental data/fitted results , 
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TABLE I: Parameters for coalescing all data points to the 
same curves in Figs. 1 and 2. 



and show _B as a function of pT in linear scale for all 
the data sets, as shown in Fig. 3 for the case of proton. 
From the figure one can see that almost all the points 
have values of B within 0.7 to 1.3, which means that the 
scaling is true within an accuracy of 30%. This is quite 
a good fit, considering the fact that the data cover about 
6 orders of magnitude. For anti-protons, the agreement 
is better than for protons. 

Now one can see that the transverse momentum dis- 
tributions for protons and anti-protons satisfy a scaling 
law. For large pt (thus large z) the scaling functions in 
Eqs. (6) and (7) behave as powers of pt, though the 
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FIG. 3: Ratio between experimental data and the fitted re- 
sults shown in Fig. 1. STAR and PHENIX data are taken 
from jSj, Symbols are the same as in Fig. 1. 
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FIG. 4: Normalized scaling distribution for protons produced 
at mid-rapidity and very forward direction in Au-|-Au col- 
lisions at RHIC with the scaling variable u. STAR and 
PHENIX data are taken from ^, ^ and BRAHMS data from 



expressions are not in powers of z ot pt- The scaling 
functions in Eqs. (6) and (7) depend on the choices of A 
and K for the case with centrality 0-12% for STAR data. 
With the variable u defined in Eq. ([3]) this dependence 
can be circumvented. (z)'s for protons and anti-protons 
are 1.14 and 1.08, respectively, with integration over z in 
the range from to 12, roughly corresponding to the pt 
range measured by STAR. The normalized scaling func- 
tions for protons and anti-protons can be obtained 
easily from Eqs. (6) and (7) and are shown in Figs. 4 
and 5, respectively together with scaled data points as 
in Figs. 1 and 2. A simple parameterization for the two 
normalized scaling functions in Figs. 4 and 5 can be given 



as follows 

= 0.064exp(13.6w-16.67w2-f 3.6?;^) , 
*p(u) = 0.086 exp(12.41w- 15. 31i;2 + 3.16w^) 

with V = ln(l + u). 
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FIG. 5: Normalized scaling distribution for anti-protons pro- 
duced at mid-rapidity and very forward direction in Au-|-Au 
collisions at RHIC with the scaling variable u. STAR and 
PHENIX data are taken from 0, i] and BRAHMS data from 

As in the case for pion distributions, one can also in- 
vestigate the Pt distributions of protons and anti-protons 
in non-central rapidity regions in Au+Au collisions. The 
only data set we can find is from BRAHMS 10] at rapid- 
ity y = 2.2 and 3.2 with centrality 0-10%. It is found that 
the BRAHMS data can also be put to the same scaling 
curves, as shown in Figs. 4 and 5. The values of corre- 
sponding parameters A and K are also given in TABLE 
I. Thus the scaling distributions found in this paper may 
be valid in both central and very forward regions for pro- 
tons and anti-protons produced in Au+Au collisions at 
RHIC at ^i]^=200 GeV. 

Now one can ask for the difference between the scaling 
functions for protons and anti-protons. After normaliza- 
tion to 1 the difference between the scaling distributions 
'J'(u) for protons and anti-protons is shown in Fig. 6. In 
log scale the difference between the two scaling functions 
is invisible at low u. To show the difference clearly a ra- 
tio r — \E'p(w)/^'p(u) is plotted in the inset of Fig. 6 as 
a function of u. The increase of r with u is in agreement 
qualitatively with data shown in where it is shown that 
p/p decreases with px monotonically. The difference in 
the two scaling functions can be understood physically. 
In Au+Au collisions there are much more quarks it, d 
than u and d in the initial state. In the central region in 
the state just before hadronization, more u and d quarks 
can be found because of the nuclear stopping effect in 
the interactions. As a consequence, more protons can 
be formed from the almost thermalized quark medium 
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than anti-protons in the small px regime. Experimen- 
tal data show that in low pt region the yield of anti- 
proton is about 80% that of protons in central Au+Au 
collisions at RHIC. This difference contributes to the net 
baryon density in the central region in Au+Au collisions 
at RHIC. On the other hand, in the large pr region, pro- 
tons and anti-protons are formed mainly from fragmen- 
tation of hard partons produced in the QCD interactions 
with large momentum transfer. As shown in [ll[, the 
gluon yield from hard processes is about five times that 
of u and d quarks. The fragmentation from a gluon to 
p and p is the same. The amount of u, d quarks from 
hard processes is about 10 times that of u, d when the 
hard parton's transverse momentum is high enough. It 
is well-known that the fragmentation function for a gluon 
to p or p is much smaller than that for a u or d (u or d) 
to p (jj) because of the dominant valence quark contribu- 
tion to the latter process. As a result, the ratio of yields 
of proton over anti-proton at large pt is even more than 
that at small px- After normalizing the distributions to 
the scaling functions the yield ratio of proton over anti- 
proton increases approximately linearly with u when u is 
large. It should be mentioned that no such difference for 
TT"*", 7r~ and tt'^, because they all are composed of a quark 
and an antiquark. 




FIG. 6: Comparison between the scaling functions for pro- 
tons and anti-protons produced at mid-rapidity in Au-|-Au 
collisions at RHIC with the scahng variable u. The inset if 
for the ratio \['p(u)/»I'p(it). 

The scaling behaviors of the px distribution functions 
for protons and anti-protons can be tested experimen- 
tally from studying the ratio of moments of the mo- 
mentum distribution, {Pt) / (pt)" = /q°° for 
n = 2, 3,4, • • •. From the determined normalized distri- 
butions, the ratio can be calculated by integrating over 
u in the range from to 12, as mentioned above, and 
the results are tabulated in TABLE II. The values of the 
ratio are independent of the parameters A and K in the 
fitting process but only on the functional form of the 
scaling distributions. If the scaling behaviors of particle 
distributions are true, such ratios should be constants in- 
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TABLE II: Ratio of moments {pt) / {pt)" for protons, anti- 
protons and pions produced in Au-|-Au collisions at RHIC. 



dependent of the colliding centralities and rapidities. For 
comparison, the corresponding values of the ratio for pi- 
ous produced in the same interactions, calculated in [6[, 
are also given in TABLE II. Because of very small dif- 
ference in the scaling distributions for protons and anti- 
protons at small u, the ratio for protons increases with 
n at about the same rate as for anti-proton for small n. 
For large n, the ratio for p becomes larger than that for p 
because of the big difference in the scaling functions for p 
and p at large u. Because of the very strong suppression 
of high transverse momentum proton production relative 
to that of pions, the ratio for pions increases with n much 
more rapidly than for p and p. 

Another important question is about the difference be- 
tween the scaling functions for protons in this paper and 
for pions in fE, 6] . Experiments at RHIC have shown that 
the ratio of proton yield over that of pion increases with 
Pt up to 1 in the region pt < 3 GeV/c and saturates 
in large px region. This behavior should be seen from 
the scaling functions for these two species of particles. 
For the purpose of comparing the scaling distributions 
we define a ratio 



(8) 



and plot the ratio i? as a function of u in Fig. 7. The ratio 
increases with u, when u is small, reaches a maximum at 
u about 1 and then decreases. Finally it decreases slowly 
to about 0.1 for very large u. The highest value of R 
is about 1.6, while the experimentally observed p over tt 
ratio is about 1 at pt '^-^ 3 GeV/c. The reason for this 
difference is two-fold. One is the normalization difference 
in defining R and the experimental ratio. Another lies in 
the different mean transverse momenta (pt)'s for pions 
and protons with which the scaling variable u is defined 
and used in getting the ratio R. 

The existence of difference in the scaling distributions 
for different species of particles produced in high energy 
collisions is not surprising, because the distributions re- 
flect the particle production dynamics which may be dif- 
ferent for different particles. In the quark recombina- 
tion models [12, 13, 14] pions are formed by combining a 
quark and an anti-quark while protons by three quarks. 
Because different numbers of (anti)quarks participate in 
forming the particles, their scaling distributions must be 
different. In this sense, our investigation results urge 
more studies on particle production mechanisms. 
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FIG. 7: Ratio between the scaling functions for 

protons and pions produced in Au+Au collisions at RHIC as 
a function of the scaling variable u. The pion scaling distri- 
bution is from 0, @] . 

IV. DISCUSSIONS 

From above investigation we have found scaling distri- 
butions for protons and anti-protons produced in Au+Au 
collisions at RHIC in both mid-rapidity and forward re- 
gion. The difference between those two scaling distri- 
butions is quite small, but they differ a lot from that for 
pions and the ratio '^^/'^tt exhibits a nontrivial behavior. 

Investigations in 0, Q and in this paper have shown 
that particle distributions can be put to the same curve 
by linear transformation on px- Though we have not 
yet a uniform picture for the particle productions in high 
energy nuclear collisions, the scaling behaviors can, in 



some sense, be cornpared to that from the string frag- 
mentation picture [15| . In that picture if there are n 
strings, they may overlap in an area of Sn and the 
average area for a string is then Sn/n. It is shown 
that the momentum distributions can be related to the 
case in pp collisions also by a linear variable change 
Pt PT{{Sn/n)AuAu/iS„/n)pp)'^/'^. Viewed from that 
picture, our fitted K gives the degree of string overlap. 
The average area for a string in most central Au+Au 
collisions is about 70 percent of that in peripheral ones 
from the values of K obtained from fitting the spectra of 
proton. If string fragmentation is really the production 
mechanism for all species of particles in the collisions, one 
can expect that the overlap degree obtained is the same 
from the changes of spectrum of any particle. In the lan- 
guage in this work, values of K are expected the same 
for pions, protons and other particles in the string frag- 
mentation picture for particle production. Our results 
show the opposite. Comparing the values of K from Q 
and this work, one can see that for pion spectrum K is 
larger for more peripheral collisions but smaller for pro- 
ton and anti-proton spectra. Our results indicate that 
other particle production mechanisms may also provide 
ways to the scaling distributions. Obviously more de- 
tailed studies, both theoretically and experimentally, are 
needed. 
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